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Technical  Approach  and  Justification 
/.  Introduction  and  Motivation 

Under  extreme  environments,  such  as  especially  low  or  especially  high  temperatures,  traditional 
liquid  lubricant  performance  is  poor.  For  example,  at  temperatures  greater  than  523  to  573  K,  oil- 
based  lubricants  typically  oxidize  or  decompose  (i).  Therefore,  to  facilitate  the  function  of  high- 
performance  mechanical  systems  that  operate  at  great  temperatures  in  naval  applications, 
identifying  robust  solid-state  lubricants  with  good  friction  and  wear  properties  is  highly  desirable. 
Inorganic  oxide,  nitride,  carbide,  and  related  ceramic  materials  represent  the  best  candidates  for 
such  solid-state  lubricants  because  of  their  high  melting  temperatures.  As  a  result  of  the  relative 
scarcity  of  work  examining  these  materials  for  tribological  performance,  there  is  a 
correspondingly  urgent  need  for  the  rapid  determination  of  the  friction  and  wear  properties  of 
candidate,  solid-state  materials. 

Data  mining  and  materials  informatics  methods  were  applied  to  the  search  for  new,  high- 
temperature  solid-state  lubricant  materials  (2).  In  particular,  it  was  used  to  generate  a  predictive 
model  that  enables  efficient  high-throughput  screening  of  inorganic  materials  with  input  from 
atomic-scale  modeling  and  experimental  testing.  This  predictive  model  was  developed  with  an 
initial  database  of  38  materials,  and  has  since  been  extended  to  approximately  500  materials. 
New  solid-state  materials  with  low  friction  coefficients  were  discovered  as  a  result.  This  is 
transformative,  because  the  materials  themselves  may  serve  as  the  basis  for  new  lubricants  at 
extreme  environments.  Additionally,  analysis  of  the  properties  of  these  low-friction  materials 
will  guide  the  future  development  of  extreme  environment  lubrication  materials. 

An  additional  discovery  that  was  identified  from  these  efforts  was  the  strong  dependence  of 
inorganic  material  wear  on  the  direction  of  sliding  and  the  quantification  of  the  activation  energy 
associated  with  the  directionality  of  wear.  This  finding  is  transformative  because  it  indicates  that 
the  performance  lifetime  of  inorganic,  solid-state  lubricants  can  be  optimized  through  material 
texture  design. 

II.  Project  Objective 

The  objective  of  this  work  was  to  develop  a  sophisticated,  materials  informatics-based  approach 
to  enable  the  rapid  and  inexpensive  identification  of  low-friction,  low-wear  inorganic  materials 
with  optimized  textures  for  use  as  solid-state  lubricants.  The  approach  entailed  a  combination  of 
data  mining  of  material  databases,  development  of  criteria  for  identification  of  lubricating 
materials,  calculations  of  potential  energy  surfaces  of  sliding  interfaces  in  candidate  materials 
using  calculations,  and  high-throughput  experimental  investigation  of  self-mated  friction. 
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III.  Application  of  Materials  Informatics  to  Identification  of  High-Temperature  Solid-State 
Lubricants 

Data  mining  and  statistical  analysis  techniques  were  used  as  a  means  of  determining  which 
materials  properties  have  the  most  effect  on  friction. 

III.  A  Material  Dataset 

An  appropriate  dataset  complete  with  a  variety  of  materials  and  a  broad  assortment  of  materials 
properties  was  required  to  enable  the  use  of  materials  informatics  in  this  task.  From  a  publication 
by  Erdemir  in  2000  (5),  it  is  shown  that  there  is  an  inverse  linear  relationship  between  a 
material’s  coefficient  of  friction  and  its  ionic  potential  which  is  defined  as  the  ratio  of  the 
material’s  cation  charge  to  the  cation  radius.  Based  on  this  understanding,  we  hypothesized  that 
other  chalcogenide  materials  would  be  similar  to  oxides  and  that  other  intrinsic  materials 
properties  would  also  display  correlations  with  friction  coefficient;  therefore,  the  material  dataset 
used  in  our  prior  work  is  comprised  of  a  variety  of  binary  chalcogenides  along  with  a  handful  of 
other  non-chalcogenides  with  data  that  has  been  provided  by  Marchman  and  Sawyer  {4). 

Table  1  summarizes  the  list  of  materials  with  experimentally  determined  friction  coefficients 
along  with  each  of  the  properties  that  are  included  in  the  material  dataset  in  addition  to 
coefficient  of  friction.  These  material  properties  were  selected  for  the  dataset  since  the  properties 
of  each  material  need  to  either  be  easily  attained  or  easily  calculated.  If  the  properties  within  the 
empirical  formula  are  too  difficult  to  acquire  or  determine,  it  will  be  exceedingly  difficult  to 
predict  the  friction  coefficients  for  a  database  of  potentially  thousands  of  materials.  The 
descriptions  of  some  of  these  properties  are: 

•  Percent  ionicity  is  a  measure  of  the  amount  of  ionic  character  of  the  bonds  within  the 
crystal 

•  Madelung  constant  is  a  property  that  is  unique  to  a  particular  crystal  structure  and  is  used 
as  a  means  of  calculating  the  electrostatic  potential  energy 

•  Electrostatic  potential  energy  is  a  measure  of  the  interaction  between  point  charges  within 
the  crystal  structure 

•  Interplanar  spacing  is  the  unidirectional  distance  between  planes  of  atoms  where  cleavage 
is  most  likely  to  occur 

•  Rij  value  is  the  cation-anion  bond  length  within  the  crystal 

•  Electronegativity  difference  is  taken  to  be  the  difference  between  the  electronegativities 
of  a  single  anion  and  a  single  cation  as  is  seen  in  the  formula  for  calculating  percent 
ionicity 
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Table  1  Material  dataset  with  16  properties  and  38  materials  used  to  develop  predictive  model  for  friction  coefficient. 


Chemical 

formula 

Structure/  phase 

Coefficient 
of  friction® 

Mohs 

hardness*’ 

Formal 

cation 

charge 

Cation 

radius 

(A)' 

Ionic 

potential 

Percent 

ionicity 

(%) 

Madelung 

constant 

Electrostatic 

potential 

(eV/atom) 

MgO 

Periclase 

0.425  (5) 

5.50 

2 

0.72 

2.778 

67.833 

-1.747 

-23.904 

SiOa 

Quartz 

0.449 

7.00 

4 

0.40 

10.000 

44.728 

-1.474 

-52.773 

AI2O3 

Corundum 

0.400  (5) 

9.00 

3 

0.54 

5.556 

56.709 

-1.216 

-28.334 

ZnO 

Zincite 

0.700  (3) 

4.00 

2 

0.74 

2.703 

55.113 

-1.642 

-23.885 

CuO 

Tenorite 

0.400  (5) 

3.50  (6) 

2 

0.77 

2.597 

44.728 

-1.365 

-20.199 

FeO 

Wustite 

0.600  (3) 

5.00  (6) 

2 

0.55 

3.636 

47.692 

-1.747 

-23.350 

M0O3 

Molybdite 

0.235  (J) 

3.50  (7) 

6 

0.69 

8.696 

33.608 

-1.392 

-61.521 

NiO 

Bunsenite 

0.500  (3) 

5.50 

2 

0.69 

2.899 

44.302 

-1.747 

-24.150 

V2O5 

Shcherbinaite 

0.310(3) 

3.25  (S) 

5 

0.79 

6.329 

55.914 

-1.486 

-58.475 

Ti02 

Rutile 

0.450  (3) 

6.20 

4 

0.86 

4.651 

59.445 

-1.600 

-47.076 

Sn02 

Cassiterite 

0.500  (3) 

6.50 

4 

0.69 

5.797 

42.166 

-1.600 

-44.890 

Zr02 

Baddeleyite 

0.500  (5) 

6.50 

4 

0.72 

5.556 

67.144 

-1.660 

-43.753 

Ag2S 

Acanthite 

0.101 

2.30 

1 

1.15 

0.870 

10.024 

-1.576 

-8.921 

WS2 

Tungstenite 

0.043  (P) 

2.50  (70) 

4 

0.60 

6.667 

1.203 

-1.283 

-30.666 

PbS 

Galena 

0.202 

2.50 

2 

1.19 

1.681 

1.550 

-1.747 

-16.957 

CU2S 

Chalcocite 

0.315 

2.80 

1 

0.77 

1.299 

10.917 

-1.567 

-9.791 

M0S2 

Molybdenite 

0.220 

1.30 

4 

0.69 

5.797 

4.314 

-1.283 

-30.486 

FeS2 

Pyrite 

0.200 

6.30 

2 

0.55 

3.636 

13.118 

-0.791 

-10.070 

ZnS 

Sphalerite 

0.527 

3.80 

2 

0.74 

2.703 

19.445 

-1.637 

-20.141 

®  All  friction  coefficients  from  tribometry  experiments  performed  by  Marchman  and  Sawyer  (4)  unless  otherwise  stated 
All  Mohs  hardness  values  from  CRC  Handbook  {11)  unless  otherwise  stated 
All  ionic  radii  from  Gersten  and  Smith  (72) 
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Table  1  Continued 


Chemical 

formula 

Structure/  phase 

Coefficient 
of  friction® 

Mohs 

hardness'^ 

Formal 

cation 

charge 

Cation 

radius 

(A)' 

Ionic 

potential 

Percent 

ionicity 

(%) 

Madelung 

constant 

Electrostatic 

potential 

(eV/atom) 

SbaSs 

Stibnite 

0.300 

2.00 

3 

0.76 

3.947 

6.782 

-1.551 

-25.842 

CdS 

Greenockite 

0.370 

3.30 

2 

0.95 

2.105 

17.965 

-1.642 

-18.681 

NiS 

Millerite 

0.240 

3.30 

2 

0.69 

2.899 

10.616 

-1.626 

-20.321 

MoSe2 

Drysdallite 

0.060  (75) 

2.00  {14) 

4 

0.69 

5.797 

3.731 

-1.283 

-29.257 

ZnSe 

Stilleite 

0.490 

5.00  (70) 

2 

0.74 

2.703 

18.331 

-1.637 

-19.222 

GaSe 

Pbs/mmc 

0.230  (75) 

2.00  (76) 

2 

0.62 

3.226 

12.794 

-1.039 

-12.054 

CoSe 

Freboldite 

0.280 

2.75  (77) 

2 

0.65 

3.077 

10.616 

-1.706 

-19.832 

CuaSe 

Berzelianite 

0.490 

2.70  (70) 

1 

0.77 

1.299 

10.024 

-1.554 

-8.855 

PbSe 

Clausthalite 

0.190 

2.75  (70) 

2 

1.19 

1.681 

1.203 

-1.747 

-16.375 

CdTe 

Zinc  Blende 

0.718 

3.00  (76) 

2 

0.95 

2.105 

4.115 

-1.637 

-16.810 

Nile 

Imgreite 

0.280 

4.00  {18) 

2 

0.69 

2.899 

0.898 

-1.706 

-18.566 

GaAs 

Zinc  Blende 

0.405 

4.50  (76) 

3 

0.62 

4.839 

3.365 

-2.455 

-43.357 

CaF2 

Fluorite 

0.372 

4.00 

2 

1.00 

2.000 

89.140 

-0.839 

-10.224 

BaF2 

Frankdicksonite 

0.392 

2.50  (6) 

2 

1.35 

1.481 

90.810 

-0.839 

-9.009 

MgF2 

Sellaite 

0.429 

5.00 

2 

0.72 

2.778 

83.174 

-0.801 

-11.583 

NaCl 

Halite 

0.303 

2.00 

1 

1.02 

0.980 

71.155 

-0.873 

-4.461 

KCl 

Sylvite 

0.319 

2.00 

1 

1.38 

0.725 

74.561 

-0.873 

-3.999 

KBr 

Rock  Salt 

0.379 

1.50  (76) 

1 

1.38 

0.725 

68.174 

-0.873 

-3.813 

YPO4 

Xenotime 

0.357 

4.50 

3“ 

0.90“ 

3.333“ 

52.884 

-2.804“ 

-51.689 

^  All  friction  coefficients  from  tribometry  experiments  performed  by  Marchman  and  Sawyer(¥)  unless  otherwise  stated 
*’  All  Mohs  hardness  values  from  CRC  Handbook  {11)  unless  otherwise  stated 
®  All  ionic  radii  from  Gersten  and  Smith  (72) 

**  Value  specific  to  yttrium  ion 
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Table  1  Continued 


Chemical 

formula 

Interplanar 
spacing  (A) 

Rij  distance 

(A) 

Melting 

temperature 

(K^ 

ENof 

cation^ 

ENof 

anion^ 

EN 

difference^ 

Density 

(g/cc)® 

Molar 

weight 

(g/mol)® 

MgO 

2.106 

2.106 

3098 

1.31 

3.44 

2.13 

3.600 

40.304 

Si02 

1.500 

1.610 

1995 

1.90 

3.44 

1.54 

2.648 

60.084 

AI2O3 

1.327 

1.855 

2327 

1.61 

3.44 

1.83 

3.990 

101.961 

ZnO 

1.796 

1.981 

2247 

1.65 

3.44 

1.79 

5.600 

81.380 

CuO 

1.277 

1.948 

1500 

1.90 

3.44 

1.54 

6.310 

79.545 

FeO 

2.155 

2.155 

1650 

1.83 

3.44 

1.61 

6.000 

71.844 

M0O3 

2.102 

1.956 

1075 

2.16 

3.44 

1.28 

4.700 

143.960 

NiO 

2.084 

2.084 

2230 

1.91 

3.44 

1.53 

6.720 

74.693 

V2O5 

2.303 

1.831 

954 

1.63 

3.44 

1.81 

3.350 

181.880 

Ti02 

1.983 

1.958 

2116 

1.54 

3.44 

1.90 

4.170 

79.866 

Sn02 

2.057 

2.054 

1903 

1.96 

3.44 

1.48 

6.850 

150.709 

Zr02 

1.290 

2.187 

2983 

1.33 

3.44 

2.11 

5.680 

123.223 

Ag2S 

2.072 

2.546 

1098 

1.93 

2.58 

0.65 

7.230 

247.801 

WS2 

3.124 

2.411 

1523 

2.36 

2.58 

0.22 

7.600 

247.970 

PbS 

2.968 

2.968 

1386 

2.33 

2.58 

0.25 

7.600 

239.300 

CU2S 

1.427 

2.306 

1402 

1.90 

2.58 

0.68 

5.600 

159.157 

M0S2 

2.980 

2.425 

1458  {19) 

2.16 

2.58 

0.42 

5.060 

160.090 

FeS2 

1.464 

2.264 

1444  {19) 

1.83 

2.58 

0.75 

5.020 

119.975 

ZnS 

1.913 

2.342 

1973 

1.65 

2.58 

0.93 

4.040 

97.440 

®  All  melting  temperature  and  density  values  from  CRC  Handbook  {20)  unless  otherwise  stated 
^  Electronegativity  values  on  the  Pauling  scale 

®  Molar  weights  according  to  the  National  Institute  of  Standards  and  Technology 
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Table  1  Continued 


Chemical 

formula 

Interplanar 
spacing  (A) 

Rij  distance 

(A) 

Melting 

temperature 

(K)* 

ENof 

cation^ 

ENof 

anion^ 

EN 

difference^ 

Density 

(g/cc)® 

Molar 

weight 

(g/mol)^ 

Sb2S3 

1.915 

2.594 

823 

2.05 

2.58 

0.53 

4.562 

339.715 

CdS 

2.599 

2.532 

1753 

1.69 

2.58 

0.89 

4.826 

144.476 

NiS 

1.635 

2.306 

1249 

1.91 

2.58 

0.67 

5.500 

90.758 

MoSe2 

3.118 

2.527 

1473 

2.16 

2.55 

0.39 

6.900 

253.880 

ZnSe 

2.004 

2.454 

1790 {21) 

1.65 

2.55 

0.90 

5.650 

144.340 

GaSe 

3.184 

2.484 

1233 

1.81 

2.55 

0.74 

5.030 

148.680 

CoSe 

1.325 

2.479 

1328 

1.88 

2.55 

0.67 

7.650 

137.890 

Cu2Se 

1.460 

2.529 

1386 

1.90 

2.55 

0.65 

6.840 

206.050 

PbSe 

3.074 

3.074 

1351 

2.33 

2.55 

0.22 

8.100 

286.200 

CdTe 

2.291 

2.806 

1365 {21) 

1.69 

2.10 

0.41 

6.200 

240.010 

Nile 

1.339 

2.648 

1133  (22) 

1.91 

2.10 

0.19 

8.384(25) 

186.293 

GaAs 

1.999 

2.448 

1511 

1.81 

2.18 

0.37 

5.318 

144.645 

CaF2 

1.366 

2.366 

1691 

1.00 

3.98 

2.98 

3.180 

78.075 

BaF2 

1.550 

2.685 

1641 

0.89 

3.98 

3.09 

4.893 

175.324 

MgF2 

1.981 

1.992 

1536 

1.31 

3.98 

2.67 

3.148 

62.302 

NaCl 

2.820 

2.820 

1073.7 

0.93 

3.16 

2.23 

2.170 

58.443 

KCl 

3.146 

3.146 

1044 

0.82 

3.16 

2.34 

1.988 

74.551 

KBr 

3.300 

3.300 

1007 

0.82 

2.96 

2.14 

2.740 

119.002 

YPO4 

2.243 

2.345** 

2268  {24) 

1.71 

3.44 

1.74 

4.800(77) 

183.877 

**  Value  specific  to  yttrium  ion 

®  All  melting  temperature  and  density  values  from  CRC  Handbook  {20)  unless  otherwise  stated 
^  Electronegativity  values  on  the  Pauling  scale 

®  Molar  weights  according  to  the  National  Institute  of  Standards  and  Technology 
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Although  most  of  the  calculated  properties  included  in  Table  1  require  only  a  simple  equation  to 
determine,  the  electrostatic  potential  energy  and  Madelung  constant  necessitates  a  different 
approach.  In  particular,  the  electrostatic  potential  energy  is  calculated  for  each  material  using  the 
Generalized  Utility  Lattice  Program  (GULP)  (25,  26)  to  determine  the  single-point  electrostatic 
lattice  energy  for  a  periodic  unit  cell  of  a  given  material,  which  is  subsequently  used  to  calculate 
the  Madelung  constant  for  the  material. 

7/7.5  Development  of  Model  for  Estimating  Coefficient  of  Friction 

We  performed  various  statistical  analysis  methods  in  order  to  develop  a  data-driven  model  for 
predicting  friction  coefficient  (27).  In  order  to  characterize  the  material  dataset  shown  in  Table  1, 
we  began  our  data-mining  approach  with  the  informatics  method  principle  component  analysis 
(PCA).  PCA  is  a  mathematical  procedure  that  allows  for  the  transformation  of  a  data  set  by 
normalizing  the  data  and  reducing  it  to  descriptive  vectors  which  are  components  that  describe 
the  variability  in  the  data.  This  transformation  is  defined  in  such  a  way  that  the  first  principle 
component  (PC)  contains  the  most  variability  in  the  data  and  the  second  PC  has  the  second  most 
variability  and  so  on.  We  were  able  to  remove  PC’s  of  lesser  significance  in  order  to  reduce  the 
dimensionality  of  the  data. 
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Fig.  1  Principal  component  analysis  a)  scores  plot  and  b)  loadings  plot  for  PCI  versus  PC2. 
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The  results  of  the  PCA  performed  on  our  material  dataset  are  given  in  Figure  1.  From  the  scores 
plot  shown  in  Figure  la,  it  can  be  seen  that  noticeable  grouping  based  on  material  chemistry 
provides  clear  distinction  between  the  oxides,  chalcogenides,  and  halides  within  our  dataset.  This 
is  an  interesting  result  since  many  of  the  properties  such  as  hardness,  melting  temperature,  molar 
weight,  and  interplanar  spacing  among  others  were  not  selected  to  reflect  the  material’s 
chemistry;  however,  the  PCA  is  still  able  to  indicate  the  separation  between  these  chemically 
similar  groupings.  Figure  lb  illustrates  the  loadings  plot  from  the  PCA.  In  this  graph,  it  is 
possible  to  evaluate  the  relationships  between  properties  as  well  as  their  respective  contribution 
to  friction. 

In  order  to  determine  which  of  the  1 5  properties  included  in  our  material  dataset  are  the  most 
important,  we  performed  a  partial  least  squares  (PLS)  regression  on  the  PCA  discussed  above. 
Specifically,  PLS  regression  is  a  type  of  linear  regression  which  takes  place  within  the 
eigenspace  of  the  PCA.  From  this  PLS  regression,  we  were  able  to  determine  that  10  of  the  15 
properties  from  Table  1  provide  the  most  importance  to  the  material  dataset.  A  property’s 
importance  is  defined  as  the  ratio  between  the  variance  in  the  data  for  that  property  and  the 
variance  in  the  data  for  all  15  properties.  These  10  important  material  properties  were  found  to 
be  Mohs  hardness,  percent  ionicity,  electrostatic  potential,  interplanar  spacing,  melting 
temperature,  cation  electronegativity,  anion  electronegativity,  electronegativity  difference, 
density,  and  molar  weight.  Since  friction  coefficient  is  a  complex  property  which  varies 
depending  on  conditions,  it  is  not  defined  as  a  finite,  absolute  value;  thus,  a  linear  regression 
model  such  as  PLS  is  not  the  optimal  method  for  developing  a  friction  model  with  the  highest 
accuracy  for  our  material  dataset.  Instead,  we  have  investigated  this  data  mining  problem 
through  a  recursive  partitioning  analysis.  With  recursive  partitioning,  rather  than  defining 
friction  coefficient  of  a  material  as  a  set  value  based  on  its  properties,  a  set  of  if-then  rules  are 
generated  which  effectively  stratify  the  materials  where  each  defined  branch  meets  a  particular 
set  of  criteria.  Because  each  branch  obtained  through  recursive  partitioning  represents  a  range  of 
specific  materials  properties,  this  type  of  model  appropriately  captures  the  variation  in 
contributions  to  a  material’s  friction  coefficient. 
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Fig.  2  Dendrogram  for  estimation  of  friction  coefficient  from  recursive  partitioning. 

Using  the  material  dataset  provided  in  Table  1,  we  performed  a  recursive  partitioning  analysis 
which  indicated  that  7  of  the  original  15  materials  properties  were  statistically  important  for  our 
material  dataset.  The  resulting  if-then  rules  from  the  recursive  partitioning  analysis  are  illustrated 
in  Figure  2.  After  applying  these  criteria  to  the  materials  within  the  dataset,  the  predicted  values 
for  the  materials’  friction  coefficients  are  shown  in  Figure  3.  The  graph  shown  is  a  plot  of  the 
predicted  friction  coefficient  versus  the  average  experimental  friction  coefficient  which  allows 
for  a  visual  estimate  of  the  model  accuracy  since  a  perfect  prediction  will  fall  along  the  45°  blue 
line.  Also,  the  error  bars  in  this  graph  represent  the  standard  deviation  applicable  to  the 
approximate  error  for  the  criteria  within  each  branch  in  Figure  2. 

As  is  indicated  in  Figure  3,  the  friction  model  we  developed  is  highly  accurate  with  an  value 
of  0.8904;  however,  before  we  can  truly  attest  to  the  accuracy  of  any  predictive  model,  it  is 
important  to  apply  an  appropriate  means  of  validation.  We  accomplished  this  through  a  leave- 
one-out  (LOO)  cross-validation  of  this  friction  model.  Through  cross-validation,  the  dataset  is 
divided  into  different  subsets.  From  these  subsets,  different  predictive  models  are  generated  with 
one  subset  always  being  removed  from  the  analysis  for  validation  of  the  predictive  model.  After 
the  LOO  cross-validation  of  this  friction  model,  the  value  only  drops  to  0.8193  which  shows 
that  the  model  we  have  developed  is  very  robust  and  accurate  for  estimating  friction  coefficients 
of  many  additional  materials  not  included  in  the  generation  of  this  model.  Based  on  the  high 
degree  of  accuracy  along  with  the  range  of  friction  criteria  captured  by  each  branch  of  this 
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friction  model,  the  recursive  partitioning  method  is  found  to  be  the  most  accurate  and  widely 
applicable  model  for  estimating  the  coefficient  of  friction  for  a  given  material. 


Fig.  3  Predicted  versus  experimental  friction  coefficient  from  recursive  partitioning.  Error  bars 
represent  the  standard  deviation  relative  to  each  branch  of  the  dendrogram  (27). 

III.C  Tribological  Database 

In  order  to  use  data  mining  to  determine  novel  candidate  solid  lubricant  materials,  a  tribological 
database  was  constructed  using  File  Maker  Pro  which  currently  contains  over  1500  materials 
which  can  be  sorted  based  on  specific  materials  properties.  File  Maker  Pro  provides  a  great  deal 
of  flexibility  as  to  the  types  and  quantity  of  information  that  can  be  added  for  each  material 
including  documents,  pictures,  movies,  and  web  URLs.  File  Maker  Pro  also  allows  for 
manipulations  of  the  graphical  user  interface  which  enables  the  database  to  be  designed  to  meet 
specific  needs  as  well  as  capabilities  to  sort  by  materials  properties.  By  combining  the 
information  within  the  File  Maker  Pro  tribological  database  with  our  predictive,  materials 
informatics  model  for  friction,  we  were  able  to  filter  the  materials  to  indicate  those  which  are 
predicted  to  exhibit  specified  coefficients  of  friction. 

III.D  Rapid  Determination  of  Friction  Coefficients  via  Macro-scale  Tribometry 

The  tribological  properties  of  19  different  interesting  minerals  were  examined  with  the  friction 
and  wear  properties  of  each  of  these  minerals  provided  in  Table  2.  A  pin-on-disk  tribometer  is 
used  to  conduct  the  tribological  test  that  provides  the  mineral’s  friction  coefficient  data.  Once  the 
test  is  complete,  a  surface  profiler  measures  the  wear  scar  depth  and  width  for  wear  rate 
calculation.  A  scanning  white  light  interferometer  is  used  to  measure  the  wear  scar  for  any 
minerals  which  the  surface  profiler  is  unable  to  track.  As  can  be  seen  from  the  data  in  Table  2, 
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the  minerals  tested  span  a  wide  range  of  friction  coefficients  and  wear  rates.  The  data  collected 
through  these  experiments  was  utilized  with  the  data  mining  methods  discussed  in  Section  1 . 

Table  2:  Friction  coefficients  and  wear  rates  for  different  minerals  determined  from  pin-on-disk 
tribometer. 


Material  Name 

Chemical 

Formula 

Average  Friction 
Coefficient 

Lowest  Friction 
Coefficient 

Highest  Friction 
Coefficient 

Standard 
Deviation  of 
Friction 

Wear  Rate 
mm^/Nm 

Acanthite 

Ag2S 

0.1010 

0.0704 

0.1431 

0.0545 

7.242E-06 

Chalcocite 

CU2S 

0.3151 

0.2470 

0.3951 

•  0.0842 

1.015E-06 

Galena 

PbS 

0.2023 

0.1699 

0.2368 

0.0422 

2.210E-05 

Malachite 

CU2C03(0H)2 

0.6608 

0.2122 

0.8241 

0.2442 

1.389E-03 

Molybdenite 

M0S2 

0.2199 

0.1624 

0.2775 

0.0935 

1.490E-04 

Pyrite 

FeS2 

0.1999 

0.1495 

0.2652 

0.0508 

1.265E-09 

Pyrophyllite 

AI2  S  ^0 -j  0  0  H2 

0.3290 

0.2322 

0.4009 

0.1487 

3.177E-02 

Realgar 

AS4S4 

0.6441 

0.1510 

0.7481 

0.1731 

3.956E-03 

Sphalerite 

ZnS 

0.5266 

0.2853 

0.6700 

0.1456 

3.034E-05 

Stilleite 

ZnSe 

0.4901 

0.1599 

0.6015 

0.1244 

1.561E-06 

Xenotime 

YPO4 

0.3565 

0.2115 

0.5678 

0.1339 

7.365E-09 

Fluorite 

CaF  2 

0.1824 

0.1229 

0.2417 

0.0447 

2.426E-08 

Sellaite 

MgFz 

0.4289 

0.1417 

0.6515 

0.1239 

2.973E-07 

Frankdicksonite 

BaF2 

0.3920 

0.2150 

0.5257 

0.1194 

6.565E-07 

Potassium 

Bromide 

KBr 

0.3787 

0.2499 

0.4294 

0.0550 

8.345E-05 

Halite 

NaCI 

0.3029 

0.2656 

0.3305 

0.0435 

2.392E-05 

Sylvite 

KCI 

0.3187 

0.2618 

0.3733 

0.0488 

9.464E-05 
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IIIE  Effect  of  Texture  on  Wear  of  Inorganic  Materials 

In  a  simple,  conceptual  model,  we  considered  the  energy  required  to  remove  an  atom  from 
within  the  surface  plane  to  sitting  on  top  of  the  plane,  neglecting  the  contribution  from  the 
migration  energy  (Fig.  4).  For  the  case  of  simple  ionic  solids,  this  energy  is  simply  the  sum  of 
the  Coulombic  potential  between  the  atom  of  interest  and  every  other  atom  within  the  surface 
plane.  This  was  performed  for  a  multiple  materials  of  different  structures  including  MgO,  NaCl, 
KBr,  KCl,  BaF2,  CaF2,  MgF2,  ZnS  and  FeS2. 


initial  state:  atom  within  surface  plane 


worn  state:  atom  on  top  of  surface  plane  and  far  from  initial  location 


neglect  the  contribution  from  the  migration  energy 


difference  in  final  and  initial  state  Is  simply  a  plane  of  atoms 


Fig.  4  Atomic  wear  of  an  ionic  solid  as  defined  as  moving  surface  atom  to  adatom  position 

The  resultant  of  this  conceptual  model  is  a  “wear  energy”  dependent  on  crystalline  structure, 
including  the  ionic  charge  of  constituents  as  well  as  lattice  spacing  and  crystalline  structure.  The 
model  is  also  dependent  on  which  crystalline  plane  is  expressed  at  the  surface. 

Comparison  of  this  wear  energy  from  the  model  with  macroscopic  wear  rates  of  ionic  solids  is 
shown  in  Figure  5.  In  general,  as  the  calculated  activation  energy  decreases,  the  wear  rate 
increases.  This  trend  includes  the  primary  factors  for  wear  of  ionic  solids:  crystalline  structure, 
elemental  charge  and  lattice  spacing.  The  results  suggest  that  future  models  should  consider  the 
removal  of  charge  balanced  groups  of  atoms,  as  this  is  more  likely. 
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Fig.  5  Measured  wear  rate  vs  activation  energy  to  “wear”  a  single  cation  as  described  in  the 
model.  All  samples  are  single  crystalline  optical  windows 

Wear  rates  of  some  single  crystalline  ionic  solids  studied  have  produced  extremely  low  wear 
rates,  rivaling  some  of  the  best  solid  lubricant  materials  available.  Detailed  analysis  of  the  wear 
of  ionic  solids  revealed  that  there  was  a  dependency  between  sliding  direction  and  wear  rate.  To 
evaluate  this  observation,  a  tribometer  was  built  under  a  scanning  white  light  interferometer  with 
a  high  precision  rotary  stage  (Fig.  6A).  The  interferometer  can  intermittently  measure  surface 
profiles  to  analyze  wear  as  a  function  of  angle  of  the  disk.  A  pin-on-disk  experiment  will  allow  a 
single  experiment  to  sweep  through  all  sliding  directions  on  a  crystalline  surface  as  well  as 
resolve  the  wear  in  that  sliding  direction. 

Wear  experiments  on  the  (001)  surface  of  rock-salts,  including  NaCl  and  MgO  (Fig.  6  C  and  D), 
revealed  that  the  material  wear  rates  have  significant  dependence  on  crystallographic  wear 
direction.  The  materials  experienced  maximum  wear  when  sliding  in  the  <100>  family  of 
directions  and  minimum  wear  when  sliding  in  the  <1 10>.  For  MgO  the  wear  rate  in  the  <100> 
direction  was  approximately  three  times  that  in  the  <11 0>  direction.  Wear  experiments  were 
performed  at  angles  relative  to  the  [100]  direction  in  six  degree  increments  revealing  wear  as  a 
sinusoidal  function  of  direction  with  90°  periodicity.  Conceptually,  sliding  in  the  100  direction 
will  readily  produce  worn  surfaces  of  lower  energy  than  sliding  in  the  110  direction.  A  negative 


14 


control  experiment  on  glass  revealed  that  these  periodicities  are  not  associated  with  errors  with 
the  tribometer  (Fig  6  B). 
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Fig.  6  Detailed  characterization  of  orientation  dependent  tribological  properties.  A)  Tribometer 
used  to  measure  wear  rate  and  friction  coefficient  as  a  function  of  sliding  direction.  B)  Wear  rate 
as  a  function  of  sliding  direction  is  shown  for  glass.  Wear  rate  as  a  function  of  sliding  direction 
on  (001)  surfaces  of  single  crystalline  c)  NaCl  and  d)  MgO.  Countersample  (pin)  was  a  3.175 
mm  diameter  alumina  ball.  Normal  force  ~  IN  and  velocity  ~  12  mm/s. 

Results  for  other  ionic  solids  suggest  a  periodicity  exists  in  more  complicated  structures  than  the 
simple  rock  salt  structure.  Practical  implications  of  this  result  suggest  that  the  crystalline  surface 
and  sliding  direction  of  a  wear  component  is  extremely  important  for  single  crystalline  ionic 
solids  and  possibly  metals.  Through  careful  material  synthesis,  alignment  and  design,  low  wear 
components  can  be  developed  for  precision  sliding  applications. 


IV.  Prediction  of  New  Solid  State  Lubricants 


The  recixrsive  portioning  model  illustrated  in  Fig.  3  has  been  applied  to  about  500  compounds 
from  the  FileMakerPro  database.  According  to  the  original  regression  tree,  45  of  these  are 
classified  as  being  very  low  friction  materials,  as  illustrated  in  Fig.  7,  with  friction  coefficients  of 
about  0.068±0.030.  Of  these,  seven  are  previously  known  to  be  solid-state  lubricants. 
Additionally,  three  of  the  materials  are  known  to  have  friction  coefficients  greater  than  the 
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predicted  values  (Ag2S=0.4,  Aul=0.25,  and  MoSe2=0.17).  This  is  thought  to  be  because  the 
regression  tree  model  assumes  single  crystal  use,  while  the  experimental  data  for  these  three 
materials  were  carried  out  on  powder  or  thin-film  samples  that  are  expected  to  have  different 
behavior.  Most  importantly,  35  have  never,  to  our  knowledge,  been  indicated  as  possible  solid- 
state  lubricants. 
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Fig.  7.  Predicted  friction  coefficients  for  materials  from  the  dataset  with  the  regression  tree 
model. 

These  newly  discovered  predicted  low-friction  materials  include  Ag2S,  BhOa,  BiFa,  MoSea,  PbO 
(litharge),  PbO  (massicot),  WSa,  AgaO,  AgaTe,  Aul,  CuTe,  ErFa,  HfF4,  HgFa,  HgO,  IrTea,  LuFa, 
NbO,  NdaOa,,  NiTea,  PbFa,  PdO,  PraOa,  PtSea,  PuF4,  Pula,  ReOa,  ReOaScBra,  ThOa,  TlCl,  TIF, 
TlFa,  Til,  UB4,  UCa,  UFa,  W2B5,  WBre,  WNa,  WTea,  and  YbFa. 

Except  for  the  parameters,  Madelung  constant  and  density,  which  are  used  in  the  regression 
model,  there  are  still  several  parameters  relevant  to  those  materials  with  the  lowest  predicted 
friction  coefficients,  as  indicated  in  Table  3.  For  example,  a  relationship  between  ionic  potential 
and  friction  coefficient  has  been  proposed  previously,  as  discussed  in  the  previous  sections.  In 
particular,  a  higher  ionic  potential  is  correlated  with  lower  friction  coefficient.  The  ionic 
potential  is  defined  as  $  =  Z/r,  where  Z  is  the  formal  charge  on  the  cation  and  r  is  the  radius  of 
the  cation.  Generally,  the  higher  the  formal  charge  or  smaller  the  radius  of  the  cation,  the  more 
the  cation  will  be  surrounded  by  anions  and  the  less  it  will  interact  with  a  neighboring  cation. 
Based  on  this  assumption,  the  large  cationic  charge  of  mineral  compovmds  indicates  the  number 
of  anions  tends  to  be  larger  than  the  number  of  cations.  This  corresponds  to  our  predicted  result 
that  29  out  of  the  45  predicted  low-friction  materials  have  more  anions  than  cations  in  their 
stoichiometric  formula.  Only  4  out  of  the  45  have  the  same  number  of  cations  and  anions  in  their 
stoichiometric  formula,  while  12  out  of  45  have  fewer  anions  than  cations. 


16 


On  the  other  hand,  the  smaller  ionic  radius  of  the  cation  creates  a  better  screening  influence  for 
cations  from  other  cations.  Correspondingly,  29  out  of  the  45  predicted  low-friction  materials 
have  a  radius  smaller  than  1  A.  If  we  calculate  the  ionic  potential  for  these  materials,  up  to  30 
materials  have  ionic  potentials  larger  than  2.5. 


Table  3:  Analysis  of  45  low  friction  materials  and  their  shared  physical  properties 


Lowest  Friction 

Yes 

No 

Ionic  potential  >2.5 

30 

15 

#  of  anions  >  #  of  cations 

29 

<4  :=12 

Radius  of  cation  <100 

29 

16 

Difference  in  electronegativity  of  cations  and  anions  <  2.0 

34 

(14  are  layered) 

11 

(9  are  layered) 

Melting  temperature  <1511  K 

34 

11 

Anion  is  a  halogen 

17 

28 

Layered  structure  (next  layer  has  the  same  charge  ) 

23 

22 

Lanthanide  cation  with  large  radius,  f-electron  chemistry 

28 

17 

The  formation  of  bond  types  is  highly  dependent  on  the  difference  of  electronegativity  (EN) 
between  the  two  elements  in  the  compoxmd.  For  instance,  the  formation  of  non-polar  covalent 
bonds,  polar  covalent  bonds,  and  ionic  bonds  are  associated  with  differences  in  electronegativity 
that  are  smaller  than  0.5,  between  0.5  and  2.0  and  larger  than  2.0,  respectively.  In  our  45  low 
friction  materials,  most  of  them  form  covalent  bonds  and  only  11  of  them  form  strong  ionic 
bonds  or  have  electronegativity  differences  larger  than  2.0.  Among  the  materials  with  strong 
ionic  bonds,  9  out  of  1 1  are  layered  structures.  Because  the  strong  ionic  bonds  establish  strong 
connections  between  atoms,  sliding  should  readily  occur  between  layers.  In  addition,  generally, 
lower  melting  temperatures  indicate  weaker  intercoimections  between  atoms  in  the  system.  In  the 
model  predictions,  34  out  of  45  materials  have  the  melting  points  lower  than  1511  K. 

We  also  note  that  there  are  properties  associated  with  15  of  the  45  materials  that  make  them 
unsuitable  candidates  for  solid-state  lubrication  that  were  not  included  in  the  regression  tree, 
namely  their  toxicity:  Bi203,  BiFs,  PbO  (Litharge),  PbO  (Massicot),  Hg2Cl2,  HgF2,  HgO,  TI2O, 
TlCl,  TiF,  TiFs,  Th02,  UB4,  UC2,  and  UF3.  Additionally,  16  of  the  45  materials  are  unsuitable 
because  their  melting  temperatures  are  lower  than  1000  K:  Aul,  BiF3,  CuTe,  Hg2Cl2,  HgF2,  HgO, 
IrTe2,  PbO  (Litharge),  Re03,  TI2O,  TlCl,  TIF,  TIF3,  TII,  WBre,  WN2.  Finally,  2  of  45  will  never 
been  used  as  solid  lubricants  due  to  their  radioactive  nature:  Bi203  and  BiF3. 

If  we  exclude  the  materials  listed  above,  the  model  predicts  24  new  candidates  for  hish 
temperature  solid  lubricant  applications  that  have  previously  been  undiscovered.  Validating 
the  model,  refining  it  when  needed,  applying  it  to  additional  compounds  in  the  database,  and 
understanding  the  rules  that  govern  the  lubricity  of  the  new  candidate  materials,  will  lead  to  the 
design  of  new  custom  lubricants  that  operate  at  extreme  temperatures.  These  team  plans  to 
undertake  these  tasks  in  the  future. 
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V.  Relevance  to  Naval  Operations 


All  modern  naval  ships  and  aircraft  rely  on  tribology  for  mobility,  weapons  and  tracking  systems, 
electronic  switching  and  control,  propulsion,  and  power  generation.  For  example,  the  turbine 
engines,  from  those  found  in  the  nuclear  reactor  driven  aircraft  carriers  and  submarines  to  more 
the  conventional  combustion  jet  turbines  on  aircraft  and  helicopters  (both  jet  and 
turbofan/turboprop),  all  rely  on  oil  delivered  additives  for  lubrication  of  the  various  gears  and 
shafts  in  motion.  Combustion  jet  engines  are  particularly  limited  in  fuel  burning  efficiency  due  to 
the  operating  temperature  thresholds  of  such  lubricants;  a  new  high  temperature  solid  lubrication 
solution  would  be  extremely  desirable. 

Tribo-corrosion,  erosion,  as  well  as  other  wear  phenomena  that  are  unique  to  naval  operations 
call  for  new  materials  that  are  impervious  to  aqueous  environments.  These  application  areas 
include  nearly  all  exposed  moving  mechanical  assemblies  on  ships,  vessels,  and  aircraft  in 
addition  to  the  various  submerged  assemblies  on  submarines  and  other  US  Navy  assets. 
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